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SUMMARY 


V' 


pThe  application  of  some  fracture  mechanics  concepts  to  the  prediction  of  the  strength 


of  pressure  vessels  in  which  cracks  have  penetrated  entirely  through  the  vessel  wall  is 
examined.  The  strength  of  cylindrical  pressure  vessels  (with  longitudinal  through- 
cracks)  is  related  to  the  geometry  of  the  cylinder,  the  crack  length,  the  material  yield 
strength,  and  fracmr^  toughness.  An  expression  for  the  crack  tip  stress  intensity  factor 
is  derived  by  considering  the  influence  on  the  stress  near  the  crack  tip  of  membrane 
stressing  in  the  cylinder  and  of  bending  caused  by  bulging  of  the  partially  constrained 
material  on  either  side  of  the  crack.  The  form  of  the  derived  expression  is  verified  and 
empirical  constants  evaluated  from  burst  data  at  -320°  and  -423°.^  of  aluminum  alloy 
(2014-T6  A1)  and  extra- low- interstitial  (ELI)  titanium  alloy  (5A1-2.  '5Sn-Ti)  tanks.  ^ 

The  analysis,  based  on  the  Irwin  fracture  theory,  provides  a  suitable  means  for  the 
correlation  of  the  strength  of  flat- sheet  specimens  with  that  of  cracked  pressurized  cylin¬ 
ders.  The  correlation  suggests  that  the  fracture  strength  of  through- cracked  cylinders 
of  any  radius,  crack  length,  and  thickness  can  be  determined  from  material  toughness, 
yield  strength,  and  a  bulge  coefficient.  Neglect  of  the  effect  of  pressure  bulging  in  the 
fracture  analysis  of  pressurized  cylinders  can  cause  a  significant  overestimate  of  the 


fracture  strength. 
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( The  presence' o^cracks  in  the  wall  of  a  pressure  vessel,  such  as  a  cryogenic  liquid 
rocket  propellant  tank,  can  severely  reduce  the  strength  of  the  structure  and  can  cause 
sudden  failure  at  nominal  tensile  stresses  less  than  the  material  yield  strength.  In  such 
cases,  cracks  form  at  inherent  flaws  in  the  structure,  extend  under  loading  to  a  critical 
size,  and  then  propagate  rapidly  through  the  material.  To  ensure  the  integrity  of  a  struc¬ 
ture  when  a  flaw  is  present  requires  that  the  designer  understand  and  adequately  consider 
the  mechanics  of  fracture,  particularly  the  relation  between  fracture  load  and  the  flaw 
and  structural  geometry.  / 


G.  R.  Irwin  has  developed  a  useful  fracture  analysis  that  relates  the  stress  distribu¬ 
tion  near  the  tip  of  a  crack  to  the  Griffith  energy  criteria  for  unstable  crack  growth 
(ref.  1).  This  analysis  relates  the  fracture  load  to  the  size  of  the  crack,  the  geometry  of 
the  structure,  and  a  single  material  fracture  toughness  parameter.  The  Irwin  analysis 
has  been  applied  to  various  structural  and  crack  configurations  including  tensioned  flat 
plates  having  through  cracks  (ref.  2)  or  surface  cracks  (ref.  3),  edge- cracked  beams  in 
bending  and  tension  (ref.  4),  and  surface- cracked  pressurized  cylinders  (ref.  5).  Al¬ 
though  a  fracture  analysis  of  a  pressurized  cylinder  in  which  a  crack  has  penetrated  en¬ 
tirely  through  the  vessel  wall  is  not  available,  there  is  evidence  that  the  strength  of 
through- cracked  pressurized  cylinders  is  less  than  that  of  a  flat  sheet  having  the  same 
crack  size.  This  behavior  was  discussed  by  Peters  and  Kuhn  in  reference  6,  which  re¬ 
ports  results  of  burst  tests  of  through- cracked  aluminum  cylinders.  .j-  < 

pfhe  present  work  examines  this  apparent  reduction  in  cylinder  burM  strength  through 
the  application  of  fracture  mechanics  concepts.  The  strength  reduction  is  considered  to 
be  due  to  an  increase  in  the  crack  tip  stress  intensity  caused  by  pressure  bulging  near 
the  crack.  The  stress  intensity  factor  associated  with  bulging  is  derived  by  a  dimensional 
analysis,  and  the  Griffith-Irwin  fracture  criterion  is  employed  to  relate  circumferential 
stress  at  fracture  to  crack  length,  cylinder  radius,  material  fracture  toughness,  yield 
strength,  and  an  empirical  bulge  coefficient.  Data  from  burst  tests  of  cylinders  of 
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2014-T6  A1  and  ELI  5A1-2.  5Sn-Ti  at  -320'°  and  -423°  F  verify  the  form  of  the  expression 


SYMBOLS 
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2a  crack  length,  in. 

2a  critical  crack  length,  in. 

2a  crack  length  indicated  by  continuity  gage,  in. 

2aQ  initial  crack  length,  in. 

C  dimensionless  bulge  coefficient 

C.  set  of  experimental  bulge  coefficient  values 
E  Young's  modulus,  psi 

fractional  error  in  fracture  toughness  value  due  to  material  inconsistencies 
K  stress  intensity  factor,  psi-^inT 

stress  intensity  factor  associated  with  pressure  bulging,  psi  -^im 
critical  stress  intensity  factor,  fracture  toughness,  psi 
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nominal  fracture  toughness  (based  on  initial  crack  length),  psi-^m 

stress  intensity  factor  associated  with  membrane  hoop  stressing  in  pressurized 
cylinder,  psi-^/im 
dimensionless  coefficient,  eq.  (2) 

pressure,  psi 
radius  of  cylinder,  in. 

coordinate  distance  from  crack  tip  in  crack  plane,  in. 

distance  from  crack  tip  to  location  of  failure  in  continuity  gage,  in. 

length  of  plastic  zone  at  crack  tip,  in. 

stress  concentration  factor  for  cracked  flat  sheet 

stress  concentration  factor  for  cracked  cylinder 

sheet  thickness,  in. 

sheet  width,  in. 

least- squares  method  weighting  factors 

dimensionless  exponents,  eq.  (2) 

failure  strain  of  continuity  gage,  in.  /in. 

strain  ahead  of  crack  and  normal  to  crack  plane,  in.  /in. 

strain  ahead  of  crack  and  normal  to  crack  plane  at  elastic-plastic  interface, 
in.  /in. 

uniform  gross  stress  acting  normal  to  plane  of  crack,  psi 

critical  uniform  gross  stress  acting  normal  to  plane  of  crack,  psi 

membrane  hoop  stress  in  cylinder,  psi 

critical  membrane  hoop  stress,  psi 

membrane  longitudinal  stress  in  cylinder,  psi 

uniaxial  yield  strength,  psi 

yield  strength  in  2:1  biaxial  stress  field,  psi 

stress  ahead  of  crack  and  normal  to  crack  plane,  psi 

surface  tensile  bending  stress  normal  to  crack  plane  ahead  of  crack,  psi 

stress  ahead  of  crack  and  normal  to  crack  plane  at  elastic- plastic  interface,  psi 
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ANALYSIS 


Fracture  Criterion 

The  stress  distribution  near  a  crack  in  a  structure  can  be  described  by  a  linear  elas¬ 
tic  stress  analysis,  provided  that  plastic  deformation  is  confined  to  a  small  region  near 
the  tip  of  the  crack.  G.  R.  Irwin  (ref.  7)  and  M.  L.  Williams  (ref.  8)  have  shown  that 
the  stress  near  the  end  of  a  crack  in  a  linear- elastic  material  is  inversely  proportional 
to  the  square  root  of  the  distance  from  the  crack  tip.  For  example,  the  stress  ahead  of 
the  crack  and  normal  to  the  crack  plane,  a  ,  varies  with  distance  r  from  crack  tip  by 
a  =  —  for  small  r.  The  constant  of  proportionality  in  this  relation  is  the  stress  in- 

^  -vfc 

tensity  factor  K,  which  expresses  the  magnitude  of  the  stress  field  near  the  crack  tip. 

The  stress  intensity  factor  is  a  function  of  the  applied  load,  the  crack  size  and  shape, 
and  the  geometry  of  the  structure. 

Irwin  (ref.  1)  has  related  the  stress  intensity  factor  to  the  rate  of  release  of  elastic 
strain  energy  with  crack  extension,  which  according  to  the  Griffith  theory  of  fracture 
reaches  a  critical  value  at  fracture.  Consequently,  a  fracture  criterion  is  formulated  in 
terms  of  the  stress  intensity  factor;  Fracture  occurs  when  K  attains  a  critical  value 
K  ,  usually  termed  the  fracture  toughness  of  the  material.  Values  of  K  can  be  ob- 
tained  from  laboratory  tests  of  suitable  precracked  specimens  by  methods  discussed  in 
detail  in  references  2  and  4.  By  use  of  the  Irwin  fracture  criterion  the  strength  of  a 
structure  containing  a  crack  can  be  determined  from  the  fracture  toughness  of  the  mate¬ 
rial  and  the  dimensions  of  the  crack  and  of  the  structure. 

Application  of  this  criterion  to  through- cracked  cylindrical  pressure  vessels  re¬ 
quires  an  expression  for  the  stress  intensity  factor  which  includes  the  special  geometri¬ 
cal  and  loading  conditions  pertinent  to  this  kind  of  structure.  If  such  an  expression  were 
known,  the  strength  of  a  cracked  cylinder  could  be  determined  from  its  dimensions,  the 
length  of  the  crack,  and  the  toughness  of  the  material. 

Stress  Intensity  Factor  for  a  Through -Cracked  Cylindrical  Pressure  Vessel 

The  stress  intensity  near  the  end  of  a  through  crack  in  a  pressurized  cylinder  is 
considered  to  be  influenced  by  two  distinct  loading  mechanisms  -  membrane  stressing  in 
the  cylinder  wall  and  bending  caused  by  bulging  of  the  partially  constrained  material  on 
either  side  of  the  crack.  The  following  analysis  examines  the  stress  intensity  associated 
with  each  of  these  loading  effects  for  the  case  of  a  longitudinal  crack  of  length  2a  lo¬ 
cated  in  a  cylinder  as  shown  in  figure  1.  The  segment  of  the  cylinder  wall  which  contains 
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the  crack  is  represented  in  figure  2(a)  as  a  sheet  subjected  to  uniform  lateral  pressure 
p  and  normal  stresses  and  or^ ,  corresponding  to  hoop  and  longitudinal  stresses,  re¬ 
spectively.  The  stress  intensity  factor  Kj^  due  to  hoop  stress  is  assumed  to  be  that 
of  a  centrally  cracked  wide  sheet  subjected  to  a  remote  uniform  tensile  stress.  Accord¬ 
ing  to  reference  3,  can  be  expressed  by 

Kh  =  CTj^  (1) 

The  contribution  of  the  longitudinal  stress  cr^  to  the  stress  intensity  is  considered  to  be 
negligible  as  discussed  in  reference  2. 

The  influence  of  lateral  pressure  on  the  stress  intensity  can  be  inferred  by  consider- 
ir^  the  effects  of  pressure  on  the  material  on  either  side  of  the  crack.  Pressure  against 
the  partially  constrained  portion  of  the  tank  wall  near  the  crack  would  cause  deflection 
and  rotation  of  the  crack  edge  as  shown  in  figure  2(b).  Since  the  rotation  is  zero  at  the 
crack  tip,  a  restraining  bending  moment  occurs  ahead  of  the  crack  thereby  introducing  a 
bending  stress  normal  to  the  crack  plane.  The  surface  tensile  bending  stress  ahead  of 
the  crack,  designated  (Tyj^,  is  assumed  to  be  a  function  of  distance  from  the  crack  tip  r, 
pressure  p,  sheet  thickness  t,  Young’s  modulus  E,  and  crack  length  2a.  In  the  ab¬ 
sence  of  an  explicit  relation  between  and  this  set  of  variables,  it  is  assumed  that 
Uyij  close  to  the  crack  tip  can  be  expressed  in  the  functional  form 

cryij(p,  E,  a,  t,  r)  =  kp“ E^a>' t^r^  (2) 

where  k  is  a  dimensionless  coefficient  and  a,  /3,  y,  6,  and  A  are  dimensionless  expo¬ 
nents.  Correct  dimensionality  of  equation  (2)  requires  that  /3  =  1  -  a  and  y  =  -(6  +  A) 
or 


‘^yb 


(3) 


The  exponents  a  and  A  can  be  evaluated  from  the  following  known  aspects  of  the 
mechanics  of  the  problem.  Since  the  material  is  linearly  elastic,  the  bending  stress  is 
proportional  to  the  pressure,  whereby  in  equation  (3),  a  =  1.  In  all  linear- elastic  crack 
problems  the  stress  near  the  crack  tip  varies  inversely  with  the  square  root  of  the  dis¬ 
tance  from  the  crack  tip  (refs.  7  and  8).  Thus  A  must  equal  -1/2.  For  later  conven¬ 
ience,  let  -6  equal  a  dimensionless  number  p.  Equation  (3)  can  now  be  written  as 


5 


(4) 


The  stress 
reference  2  by 


cr. 


near  the  tip  of  a  crack  is  related  to  the  stress  intensity  factor  K  in 


K  =  Oy  V2irr  (5) 

Consequently,  the  bending  stress  CTyj^  close  to  the  crack  tip  can  be  related  to  a 
stress  intensity  factor  associated  with  pressure  bulging  by 

Kb  =  C^yb  (6) 

Then  from  equations  (4)  and  (6)  this  bulge  stress  intensity  factor  is  related  to  pres¬ 
sure,  crack  length,  and  sheet  thickness  by 


=  Cp  (7) 

where  C  is  a  dimensionless  coefficient  equal  to  k  V2. 

For  cases  in  which  both  membrane  stressing  and  pressure  bulging  occur  near  a 
crack,  the  total  stress  intensity  is  the  sum  of  the  stress  intensities  for  each  of  the  load¬ 
ing  conditions,  since  stress  intensity  factors  can  be  superimposed  for  a  given  crack 
length.  In  the  special  case  of  a  pressurized  cylinder  having  a  crack  oriented  parallel  to 
its  axis,  the  total  stress  intensity  factor  K  is  the  sum  of  the  stress  intensity  factor  due 
to  hoop  membrane  stressing  Kj^  and  that  due  to  bulging  K^;  that  is, 

K  =  Kj^  +  (8) 


From  equations  (1),  (7),  and  (8), 


K  = 


(9) 


Equation  (9)  is  a  stress  intensity  factor  for  a  cracked  plate  element  under  uniform  tensile 
stress  and  lateral  pressure  p.  When  this  element  is  considered  to  be  a  segment  of 
the  cylinder  wall,  p  =  Cj^t/R  for  a  cylinder  of  radius  R,  and  equation  (9)  becomes 
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At  fracture,  K  equals  K^,  2a  equals  the  critical  crack  length  2a^,  and  equals  the 
critical  hoop  stress  which,  from  equation  (10),  is 


K. 


^hc  = 


(11) 


It  should  be  noted  that  as  R  becomes  very  large,  equation  (11)  reduces  to  the  expression 
for  the  critical  stress  on  a  wide  flat  sheet,  which  according  to  reference  2  is 


cr 


c 


(12) 


In  this  expression  the  crack  length  is  usually  altered  to  account  for  the  effect  of  plastic 
deformation  at  the  tip  of  the  crack  by  replacing  a^  with  a^  +  — —  (ref.  9)  where  (t„ 

C  C  o  i 

2iTffY 

is  the  uniaxial  yield  strength.  If  a  similar  plasticity  correction  is  used  in  the  corre¬ 
sponding  square  root  term  in  equation  (11),  the  critical  hoop  stress  becomes 


‘^hc  ■ 


K 


c 


(13) 


where  is  the  yield  strength  in  a  2:1  biaxial  stress  field.  This  expression  relates 

the  critical  hoop  stress  in  a  cracked  pressurized  cylinder  to  crack  length,  wall  thickness, 
cylinder  radius,  yield  strength,  and  fracture  toughness.  In  equation  (13)  the  crack  length 


in  the  bulge  term 


was  not  modified  for  plastic  deformation  because  the  re¬ 


sulting  expression  would  indicate  a  finite  bulge  effect  at  zero  crack  length. 

The  form  of  equation  (13)  can  be  verified  and  the  values  of  p  and  C  obtained  from 
experimental  data  from  burst  tests  of  precracked  cylinders.  When  equation  (13)  is 


rewritten  as 


where 


(14) 


the  parameters  C  and  p  can  be  separated  by  taking  the  logarithm  of  each  side  of  equa¬ 
tion  (14): 


& 

In  i// =  (p  -  l)ln -5  +  In  C  (15) 

t 

This  expression  predicts  a  linear  relation  between  i//  and  a^/t  on  a  log-log  plot. 

Burst  tests  of  pressurized  cylinders  are  required  to  test  this  relation  and  to  evaluate 
the  parameters  ju  and  C. 


APPARATUS  AND  PROCEDURE 

The  validity  of  the  form  of  equation  (13)  was  investigated  using  data  from  burst  tests 
of  cylinders  of  four  different  materials  having  various  thicknesses,  radii,  and  crack 
lengths,  at  room  temperature,  -320*^,  and  -423*^  F.  Fracture  data  for  2014-T6  A1  and 
ELI  5A1-2.  5Sn-Ti  were  obtained  at  the  Lewis  Research  Center  from  tests  of  6-inch- 
diameter  cylinders  such  as  that  shown  in  figure  3  having  longitudinal  notches  or  cracks 
varying  in  length  from  0. 1  to  2.  0  inches.  The  2014- A1  cylinders  were  cut  from  extruded 
tubing,  heat  treated  to  the  T6  condition,  and  machined  to  a  wall  thickness  of  0.  060  inch 
(ref.  10).  The  roots  of  the  notches  in  the  aluminum  cylinders  were  machined  to  a  radius 
of  about  0.  0002  inch.  The  aluminum  nominal  fracture  toughness  values  were  obtained 
from  tensile  specimens  that  were  cut  from  flattened  tube  stock  that  had  been  heat  treated 
with  the  cylinders.  These  specimens,  as  that  shown  in  figure  4,  were  3  inches  wide, 
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0.  060  inch  thick,  and  had  1-inch  central  notches  machined  to  approximately  0.  0002  inch 
radius  at  the  root. 

The  titanium  cylinders  were  formed  from  0.  020- inch  sheet  material  with  a  single 
longitudinal  butt  weld.  Each  cylinder  was  stress-relieved  by  heating  to  1100  F  for 
2  hours.  Cracks  were  started  by  extending  half-wall-thickness  electrical  discharge 
notches  through  the  wall  by  low-stress  pressure  cycling.  Fatiguing  was  continued  until 
the  cracks  reached  the  desired  test  length.  The  3-inch-wide  titanium  toughness  speci¬ 
mens  were  made  from  the  same  sheet  stock  as  that  used  for  the  cylinders  and  were  sub¬ 
jected  to  the  same  heat  treatment  and  crack  preparation. 

Room  temperature  data  for  two  aluminum  alloys  (2024-T3  and  7075-T6  Al)  were  ob¬ 
tained  from  results  of  notched  cylinder  burst  tests  published  in  reference  6.  These  data 
represent  cylinder  radii  of  3.  6  and  14.  4  inches,  wall  thicknesses  ranging  from  0. 006  to 
0. 025  inch,  and  crack  lengths  varying  from  0. 24  to  7. 70  inches. 

The  direction  of  testing  for  both  the  2014-T6  Al  and  ELI  5A1-2.  5Sn-Ti  toughness 
specimens  corresponded  to  the  hoop  direction  in  the  cylinders.  During  loading  of  the 
toughness  specimens  and  the  cylinders,  slow  crack  extension  was  monitored  by  electrical 
resistance  foil  element  continuity  gages  (ref.  11)  placed  at  the  crack  tip.  These  gages 
consisted  of  20  elements  connected  in  parallel  and  were  mounted  normal  to  the  crack  as 
shown  in  figure  5.  As  the  crack  grew  individual  elements  broke  causing  an  abrupt  mea¬ 
surable  resistance  change.  In  the  toughness  tests,  the  initial  crack  length  2a^,  the  criti¬ 
cal  crack  length  2a.^,  and  the  critical  applied  stress  were  measured.  Fracture  toughness 
values  were  calculated  from  these  data  by  the  expression  given  in  reference  2 


where  W  is  the  sheet  width.  Because  there  were  cases  in  which  the  critical  crack 
length  was  not  known,  an  alternate  correlation  was  based  on  the  initial  crack  length,  and  a 
nominal  fracture  toughness  was  computed  by 


Although  exact  values  of  K  for  the  2024-T3  and  7075-T6  Al  sheet  used  in  the  tests  of 
reference  6  are  not  available,  typical  toughness  values  of  90  000  psiVinch  for  the  2024-T3 
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and  53  000  psi  Vinch  for  the  7075-T6  A1  were  used  in  the  computations.  These  values 
were  estimated  from  flat  sheet  strength  data  included  in  reference  12. 

In  the  burst  tests  of  2014-T6  A1  and  ELI  5A1-2.  5Sh-Ti  cylinders  measurements  were 
made  of  the  initial  crack,  the  critical  pressure,  and  the  critical  crack  length  as  indicated 
by  the  continuity  gages.  The  gage- indicated  crack  length  was  subsequently  corrected  to 
obtain  a  more  accurate  critical  crack  length  measurement,  as  discussed  in  appendix  A. 
Table  I  summarizes  pertinent  cylinder  burst  test  data. 

Values  of  biaxial  yield  strength  Oyg  of  the  aluminum  were  calculated  from  the 
Von  Mises  yield  criterion  to  be  1.  15  Oy  for  a  2:1  stress  field.  Because  some  titanium 
alloys  exhibit  significant  increase  in  biaxial  yield  strength  relative  to  uniaxial  yield 
strength,  values  of  Oyg  of  ELI  5A1-2.  5Sn-Ti  were  obtained  from  pressure  tests  of  un¬ 
cracked  cylinders. 

Data  from  the  cracked  cylinder  burst  tests  were  used  to  verify  the  expression  for 
fracture  strength  developed  in  the  ANALYSIS  section  as  given  by  equation  (15).  The  form 
of  this  expression  was  tested  and  the  empirical  parameters  ju  and  C  evaluated  by  plot¬ 
ting  experimental  values  of  against  a/t  on  a  log- log  plot.  The  values  of  were 
computed  for  each  cracked  cylinder  burst  test  using  R,  and  a^  or  a^  for  each 

test  and  or  from  the  flat  sheet  toughness  tests. 


RESULTS  AND  DISCUSSION 

The  expression  for  the  burst  strength  of  through- cracked  pressurized  cylinders 
given  by  equation  (15)  in  the  ANALYSIS  section  indicates  a  linear  relation  between  the 
terms  x}/  and  a/t  on  a  log-log  graph.  The  plots  of  experimental  values  of  xj/  against 
a/t  in  figure  6  verify  the  predicted  linear  relation.  In  addition,  for  each  case  shown  in 
figure  6  it  appears  that  x{/  is  insensitive  to  the  ratio  a/t.  The  data  indicate  that  x}/  is 
constant  over  three  decades  of  the  ratio  a/t,  that  is,  the  slope  of  the  line,  ju  -  1,  is 
zero.  This  behavior  suggests  that  in  equation  (15)  the  empirical  constant  C  equals  xj/ 
and  can  be  computed  from  the  data  using  equation  (14) 


A 


c 


R 


a 


c 


(18) 
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for  correlation  based  on  the  critical  crack  length  or  by 


for  correlation  based  on  the  initial  crack  length. 

Weighted  averages  of  the  experimental  values  of  the  empirical  bulge  coefficient  C 
were  computed  by  the  method  discussed  in  appendix  B  and  are  shown  by  solid  lines  in  fig¬ 
ures  6  and  7.  The  value  of  C  obtained  from  the  data  of  reference  6  is  the  unweighted 
mean  of  the  data.  The  number  thus  obtained  is  approximately  the  same  as  the  coefficient 
of  a/R  reported  in  reference  6  in  which  the  ratio  of  the  stress  concentration  factor  for 
a  cracked  cylinder  to  that  for  a  cracked  flat  sheet  was  found  empirically  to  be 

=  1  +  9.25. 

Su  R 

In  that  report  the  coefficient  9.  2  was  calculated  from  a  comparison  of  the  strength  of 
cracked  flat  sheets  to  the  strength  of  pressurized  cylinders  having  the  same  crack 
lengths. 

From  examination  of  figure  7  it  appears  that  the  bulge  coefficient  C  tends  to  de¬ 
crease  with  decreasing  temperature.  At  the  present  time  it  is  not  clearly  understood 
why  C  is  temperature  or  material  dependent.  Further  investigation  is  required  to  ex¬ 
plain  this  dependency. 

The  scatter  of  the  data  in  figures  6  and  7  can  be  shown  to  be  caused  by  an  error  in 
the  half  crack  length  a  of  ±0. 02  inch  and  a  total  error  in  fracture  toughness  of  ±4  per¬ 
cent  for  the  aluminum  and  ±6  percent  for  the  titanium.  The  fractional  error  in  C  due  to 
both  of  these  errors  is  derived  in  appendix  C.  The  dashed  curves  on  either  side  of  the 
weighted  mean  values  of  C  in  figure  7  represent  the  computed  deviation  in  C  due  to  the 
estimated  errors  in  crack  length  and  fracture  toughness.  Because  these  curves  bound 
almost  all  of  the  data,  the  scatter  is  accounted  for  by  the  errors  cited  previously. 

Figures  8  to  10  are  curves  of  fracture  strength  of  pressurized  cylinders  computed  as 
a  function  of  crack  length  by  the  following  expressions  obtained  from  equations  (18)  or(19) 
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(20) 


or 


The  curves  in  figures  8  to  10  were  computed  from  equations  (20)  or  (21)  using  toughness 
values  obtained  from  flat  sheet  tests  and  the  weighted  mean  C  values  from  figure  7. 

The  dashed  curves  in  these  figures  are  the  upper  and  lower  values  of  fracture  strength 
prediction  due  to  the  experimental  error  in  C  discussed  previously.  The  excellent  cor¬ 
relation  in  these  figures  suggests  that  for  a  given  material  and  temperature  the  fracture 
strength  of  cylinders  of  any  radius,  crack  length,  or  wall  thickness  can  be  determined 
from  as  few  cylinder  burst  tests  as  are  necessary  to  obtain  a  meaningful  value  of  the 
bulge  coefficient  C. 

Neglect  of  the  effects  of  pressure  bulging  will  lead  to  an  overestimate  of  the  strength, 
as  can  be  seen  in  figures  8  to  10  by  comparing  the  fracture  strength  curves  that  include 
the  bulge  term  to  those  that  do  not.  It  is  apparent  that  there  can  be  a  large  reduction  in 
the  strength  of  cracked  cylinders  due  to  pressure  bulging,  particularly  for  larger  values 

of  C-. 

R 

CONCLUSIONS 


i  /  ■  . 
t/ 


■rj 


H  The  fracture  mechanics  analysis  of  the  effects  of  pressure  bulging  near  the  crack  on 
the  stress  intensity  facto^  as  described  in  this  report|^  provides  an  explanation  for  the 
lower  strength  of  cracked  cylinders  relative  to  flat  sheets  as  well  as  a  method  for  corre¬ 
lating  the  strength  of  cracked  pressurized  cylinders  with  that  of  flat  sheets. 

Application  of  a  fracture  analysis  of  through- cracked  pressurized  cylinders  that  does 
not  include  the  effects  of  pressure  bulging  could  result  in  an  overestimate  of  the  fracture 
strength.  ! 
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frhe  close  agreement  between  the  form  of  the  derived  expression  for  fracture  strength 
and  the  experimental  data  suggests  that  the  fracture  strength  of  cylinders  of  any  radius, 
crack  length,  and  wall  thickness  can  be  determined  from  material  toughness,  yield 


strength,  and  the  bulge  coefficient  obtained  from  tests  of  a  few  cylinders.  / 


Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 
Cleveland,  Ohio,  October  1,  1965. 
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APPENDIX  A 


CALCULATION  OF  CORRECTED  CRACK  LENGTH  FROM  CONTINUITY  GAGE  DATA 

In  the  cylinder  burst  tests  and  fracture  toughness  tests,  the  crack  lengths  indicated 
by  the  foil  continuity  gages  were  found  to  be  longer  than  the  actual  crack  length  because 
the  continuity  gage  elements  failed  slightly  ahead  of  the  actual  crack  front.  Consequently, 
a  more  accurate  measurement  of  crack  length  required  a  calibration  of  the  continuity 
gages.  The  following  theoretical  calibration  relates  the  gage- indicated  crack  length  to  a 
corrected  crack  length  in  terms  of  applied  stress,  material  yield  strength,  yield  strain, 
gage  failure  strain,  and  specimen  width.  It  is  presumed  that  the  gage  elements  failed 
ahead  of  the  crack  because  the  ductility  of  the  gage  mounting  cement  was  less  than  that  of 
'the  specimen  material.  How  this  behavior  might  affect  the  gage  readings  is  shown  in  fig¬ 
ure  11,  a  sketch  of  strain  ahead  of  a  crack.  When  a  gage  element  breaks  at  a  strain  e^g, 
it  indicates  a  crack  longer  than  the  actual  surface  crack  by  an  amount  2r£g.  If  the  gage 
indicates  a  crack  of  length  2a  ,  the  actual  crack  length  is 

o 

2a  =  2ag  -  (Al) 

The  distance  is  related  to  the  length  of  the  plastic  region  r^,  the  strain  at  the 
elastic -plastic  interface  and  the  failure  strain  of  the  continuity  gage  by  a 

simplified  approximation  of  the  expression  for  strain  ahead  of  a  crack  given  in  equations 
2--43  to  2“45  of  reference  13.  Then 


Equation  (A2)  can  be  solved  for 


by  noting  that 


Thus. 


(A3) 


The  corrected  crack  length  and  the  gage- indicated  crack  length  are  related  from  equa¬ 
tions  (Al)  and  (A3)  by 
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(A4) 


For  a  plate  of  width  W,  Tp/a  is  given  in  equation  2-28  of  reference  13  as  a  function  of 
the  ratio  of  gross  normal  stress  a  to  the  stress  at  the  elastic-plastic  interface  Oyy  by 

—  arcsin  (sin  —  sec  -  1  (A5) 

a  Tra  y  W  / 

For  an  infinitely  wide  plate,  equation  (A5)  reduces  to 


(A6) 


The  ratio  a/ag  can  be  calculated  as  a  function  of  cr/ (oJ’  O'/ Oy g  when  taken 

as  the  biaxial  yield  strength)  from  equations  (A3),  (A4),  and  (A6) 


For  a/W  «  1, 


The  continuity  gage  failure  strain  e^g  was  measured  from  a  tensile  test  of  the  gage 
mounted  on  an  aluminum  tensile  specimen.  The  measured  strain  at  which  the  continuity 
gage  failed  was  taken  as  e ^g.  In  equations  (A7)  and  (A8),  the  stress  at  the  elastic- 
plastic  interface  Uyy  was  taken  as  the  2:1  biaxial  yield  strength  Uyg  for  the  reasons 
discussed  on  page  52  of  reference  13.  The  strain  at  the  elastic -plastic  interface  Cyy  is 
1  ^ 

computed  as  ^  ~  the  same  reference. 

E 
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This  theoretical  ratio  is  plotted  in  figure  12  and  is  compared  to  experimental  mea¬ 
surements  of  for  ELI  5A1-2.  5Sn-Ti  and  2014-T6  A1  flat  sheet  specimens  at 

-320°  F.  The  agreement  is  good  enough  to  permit  the  use  of  equations  (A7)  or  (A8)  to 
correct  the  continuity  gage  readings. 
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APPENDIX  B 


CALCULATION  OF  AVERAGE  BULGE  COEFFICIENT 

A  least- square  method  is  selected  to  compute  an  average  value  of  C  from  the  set  of 
experimental  values  Since  the  observed  do  not  all  have  the  same  precision,  the 

squares  of  the  deviations  are  weighted  by  factors  (ref.  14).  The  sum  of  the  weighted 
squares  is 


n 

s  =  Y,  "^i^^^i  - 


i=l 


(Bl) 


When  minimized,  equation  (Bl)  is 


n 


-2  5)  w,(Ci  -  C) 
i=l 


(B2) 


and  at  the  minimum 

C  =  -  (B3) 

n 

i=l 

If  the  C.  values  are  weighted  such  that  the  values  of  w^  are  inversely  proportional  to 
the  square  of  the  expected  error  AC  caused  by  error  in  crack  length  and  toughness  (eq. 
(C9)  in  appendix  C), 


(B4) 


and 
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n 


where  AC|  is  given  by  equation  (C9). 
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APPENDIX  C 


ERROR  IN  BULGE  COEFFICIENT 

Because  values  of  the  bulge  coefficient  computed  from  equations  (18)  and  (19)  are 
highly  dependent  upon  the  crack  length  and  the  fracture  toughness,  an  error  in  C  is  pre¬ 
sumed  to  be  due  mainly  to  error  in  these  two  parameters.  The  magnitude  of  the  largest 
difference  in  C  expected  from  changes  in  and  a  is 
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An  error  in  measurement  of  K  due  to  an  error  Aa  is  given  by  (0K  /0a  )Aa  .  In  a 

C  C  C  G  C 

flat- sheet  toughness  test  is  related  to  a^  by 


if  the  width  effect  is  ignored.  Then, 

_ c  _ _ ^ 

0a  2a^ 


(C5) 


(C6) 


When  d^5  the  fractional  error  in  due  to  material  inconsistencies  is  included,  the 

total  error  AK  is 
c 


and 


AK  =_^  Aa„  +JK^ 


(C7) 


AK  /Aa  ' 
c  /  c 


+ 


K 


2  a. 


(C8) 


where  the  subscript  T  refers  to  a  flat-sheet  toughness  test. 
The  fractional  error  in  C  as  given  by  equation  (C4)  is  then 


(C9) 


or,  when  the  correlation  is  based  on  initial  crack  length  and  nominal  fracture  toughness, 
the  fractional  error  in  C  is 
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(CIO) 


The  error  in  C  was  computed  using  an  estimated  Aa  of  0.  02  inch  and  a  toughness 
test  crack  length  of  1.  0  inch.  The  estimated  fractional  error  in  toughness  |j|  was 
taken  as  0.  04  for  the  aluminum  and  0.  06  for  the  titanium. 
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TABLE  I.  -  CYLINDER  BURST  TEST  DATA 


Temper¬ 

ature, 

Op 

Critical 

hoop 

stress, 

psi 

Initial 

crack 

length, 

2ao, 

in. 

Critical 

crack  length 

indicated  by 

continuity  gage, 

2a  , 
g' 
in. 

Corrected 

critical 

crack 

length, 

in. 

Temper¬ 

ature, 

Op 

Critical 

hoop 

stress, 

^hc’ 

psi 

Initial 

crack 

length, 

in. 

Critical 

crack  length 
indicated  by 
continuity  gage, 

in. 

Corrected 

critical 

crack 

length, 

2a^, 

in. 

2014-T6  A1 

5A1-2 

.  5Sn-Ti 

-423 

82  200 

0.  104 

(a) 

(a) 

-423 

171  500 

0.093 

(b) 

(b) 

63  400 

.250 

160  900 

.155 

39  600 

.  500 

133  300 

.190 

32  000 

.750 

121  400 

.277 

21  000 

1.000 

114  200 

.280 

19  800 

1.250 

84  600 

.  523 

13  200 

1.750 

76  000 

.482 

11  900 

2.000 

If 

63  600 

.803 

63  000 

.  941 

-320 

71  600 

0.  113 

(a) 

(a) 

61  500 

,762 

70  600 

.150 

0.450 

0.250 

51  600 

.  985 

63  700 

.200 

.430 

.282 

40  600 

1.600 

58  500 

.250 

.430 

.310 

37  700 

1.  560 

f 

52  200 

.  300 

,  520 

.  402 

47  400 

.400 

.620 

.  512 

-320 

190  400 

0.  125 

0.300 

0.  125 

40  100 

.500 

.760 

.670 

164  900 

.241 

(a) 

(a) 

30  200 

.750 

1.002 

.952 

156  600 

.229 

.420 

.250 

23  100 

1.000 

1.280 

•  1.236 

115  500 

.440 

.780 

.616 

18  600 

1.250 

1.410 

1.380 

105  100 

.453 

.  530 

.440 

14  400 

1.750 

1.940 

1. 

914 

84  900 

.733 

.960 

.880 

11  300 

2.000 

2.120 

2. 

100 

74  600 

.765 

.860 

.790 

71  800 

946 

1 

ifin 

1  080 

Room 

64  400 

0.  114 

(a) 

(a) 

66  200 

.991 

1.260 

1. 

180 

34  000 

.500 

0.670 

0. 

586 

44  100 

1.480 

1.770 

1.720 

20  600 

1.000 

(a) 

(a) 

■ 

35  900 

1.460 

1.660 

1.620 

9  700 

2.000 

2. 

160 

2. 

130 

^Slow  crack  growth  was  not  measured. 

'^No  detectable  slow  crack  growth  occurred. 


23 


Figure  4.  “  Notched  flat  sheet  fracture  toughness  specimen. 
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Figure  5.  -  Foil  element  continuity  gage  mounted  at  notch  root. 


Experimental  values  of  parameter 


Hoop  fracture  stress,  psi 


Hoop  fracture  stress,  psi 


\ 

— 

O  Experimental  data 

- Eq.  (21) 

- Estimated  data  scatter 

- Uncorrected  for  bulge  effect 


0  .2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  2.0 

Initial  crack  length,  2aQ,  in. 

(d)  Room  temperature;  nominal  fracture  toughness  =  58  400  psi  -Vinch;  yield  strength 
OyB  ^  coefficient  C  =  6.37. 


O  Experimental  data 

■ - Eq.  (20) 

- Estimated  data  scatter 

- Un corrected  for  bulge  effect 


s 

\ 

s 

--- 

.2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  2.0 

Critical  crack  length,  2a^,  in. 

(e)  Room  temperature;  critical  fracture  toughness  =  64  500  psi  Vinch;  yield  strength 
Oyg  =  78  100  psi;  bulge  coefficient  C  =  6. 51. 

Figures.  -Concluded. 


Experimental  data 
Eq.  (21) 

Estimated  data  scatter 
Uncorrected  for  bulge  effect 


(a)  Temperature,  -423°  F;  nominal  fracture  toughness  =  87  500  psi  Vinch;  yield  strength 
Oyg  =  252  000  psi;  bulge  coefficient  C  =  1.57. 


’IT 

r 

_ 

W\n 

\ 

T' 

\ 

[:i 

O  Experimental  data 

- Eq.  (21) 

- Estimated  data  scatter 

- Uncorrected  for  bulge  effect 


.2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  2.0  2 

Initial  crack  length,  2aQ,  in. 

(b)  Temperature,  -320°  F;  nominal  fracture  toughness  =  139  000  psi  -Vinch;  yield  strength 
OyB  “  222  000  psi;  bulge  coefficient  C  =  4.30. 

q3_ — I - - - 1 - - - - - - - - - - - - - ^ - - - - - - — 


O  Experimental  data 

- Eq.  (20) 

- Estimated  data  scatter 

- Uncorrected  for  bulge  effect 


Critical  crack  length,  2a^,  in. 

(c)  Temperature,  -320°  F;  critical  fracture  toughness  K^.  =  160  000  psi  -v/inch;  yield  strength 
Oy  g  =  222  000  psi;  bulge  coefficient  C  =  4. 38. 

Figure  9.  -  Fracture  strength  of  5A1-2. 5Sn-Ti  cylinders.  Radius  of  cylinder  R  =  3. 00  inches. 


elements 


Eq.  (A8)  where 
02 


(a)  5AI-2.  5Sn-Ti. 


(a)  Foil  element  continuity  gage. 


Eq.  (A8)  where 


(b)  Strain  ahead  of  crack. 

(b)  2014-T6  AL 

Figure  11,  -  Failure  of  continuity  gage  elements  ahead  of 
crack  front. 

Figure  12.  -  Continuity  gage  calibration.  Temperature, 
-320°  F. 

NASA-Langley,  1966  E-2885 

